DTRC-PAS-90-45  Differentia!  Forms  of  the  Turbomachinery  Equation  AD-A226  183 


OTIC  FILE  COPY 


David  Taylor  Research  Center 

Bethesda,  MD  20084-5000 

DTRC-PAS-90-45  July  1990 

Propulsion  and  Auxiliary  Systems  Department 

Research  and  Development  Report 


$> 


A  Differential  Turbomachinery 
Equation  with  Viscous  Correction 


Herman  B.  Urbach 


OIK 

O.  l  t.  c  rt 

SEP  o  6  1980 


Approved  for  public  release;  distribution  is  unlimited. 


CODE  Oil  DIRECTOR  OF  TECHNOLOGY,  PLANS  AND  ASSESSMENT 
12  SHIP  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

1.  DTRC  reports,  a  formal  series,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  regardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semiformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case-by-case  basis. 


NDW-DTNSRDC  5602  51  (Rev  2 - 88 ) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


l:lsf  :  i  f -i  =WH :  >  i I  a  i  V  j 


Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OECLASSIF1CATION/OOWNGRADING  SCHEDULE 


PERFORMING  ORGANIZATION  REPORT  NUMBERS) 

DTRC-PAS-90-45 


la  NAME  OF  PERFORMING  ORGANIZATION 

David  Taylor  Research  Center 


Sc.  ADDRESS  (City.  SUM.  and  HP  Coda) 

Annapolis,  Maryland  21402-5067 


Sa.  name  OF  FUNDINQ/SPONSORINQ 
ORGANIZATION 


Sc.  ADDRESS  (City,  SUM.  and  ZIP  Coda ) 


1b.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION/AVAILABIUTY  OF  REPORT 


Approved  for  public  release;  distribution  is  unlimited. 


MONITORING  ORGANIZATION  REPORT  NUMBERS) 


8b.  OFFICE  SYMBOL 
(ItappOeaHa) 

Code  2720 


7*.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  SUM.  and  ZIP  Coda) 


lb.  OFFICE  SYMBOL 
fltapplkabla ) 


PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


WORK  UNH* 
ACCESSION  NO. 


11.  TITLE  (fnctucto  Security  C/MsJfJeatforr) 

A  Differential  Turbomachinery  Equation  with  Viscous  Correction 


12.  PERSONAL  AUTHOR<8) 

Dr.  Herman  B.  Urbach 


13a.  TYPE  OF  REPORT 

Final 


13b.  TIME  COVERED 

14.  DATE  OF  REPORT  (YEAR.  MONTH.  DAY) 

1989  1990 

FROM  TO 

1990  June  30 

IS.  PAGE  COUNT 


is.  supplementary  notation  Base(j  jjpQjj  a  thesis  submitted  in  partial  fulfillment  of  the  degree  requirements  of  the  Aerospace  Department  of  the 
University  of  Maryland,  College  Park,  Md.  Professor  Everett  Jones,  graduate  thesis  advisor. _ _ 


16.  SUBJECT  TtRMQ  (Contlnt*  on  nvw  If  niciMafy  mnd  ktonlty  by  bhck  numb*) 


C08ATI  CODES 


SUB-GROUP 


Differential  turbomachinery  equation;  Viscous  correction  for  turbomachinery  equation 


i.  ABSTRACT  (Conflhu#  on  nwfM  If  mcmary  and  ktonttty  by  bkxk  numb*) 

A  differential  turbomachinery  equation  describing  the  energy  transfer  between  a  Quid  and  any  body  moving  in  that  fluid  was 
derived.  The  derivation  is  based  upon  the  Coriolis  form  of  the  Navier-Stokes  Equations.  A  differential  equation  for  the  total  relative  ro- 
thalpy  was  also  obtained.  The  equalities  contain  a  rigorous  viscous  correction  for  the  total  enthalpy  and  rothalpy. 

The  differential  equation  defines  the  substantial  derivative  of  the  total  enthalpy  at  any  point  in  a  fluid.  This  substantial  derivative 
represents  the  energy  transfer  rate  into  or  out  of  the  fluid  at  the  point  In  order  to  obtain  the  total  enthalpy  transfer  between  the  fluid  and 
the  body,  it  is  necessary  to  integrate  over  the  entire  domain.  For  ideal  flow  regimes,  the  integration  domain  may  be  restricted  entirely 
within  the  rotor.  For  real  fluids,  viscous  coupling  requires  that  regions  outside  the  blade  area  be  considered. 

On  integration  of  the  differential  equations,  a  form  of  the  Euler  Turbomachinery  Equation  with  viscous  correction  is  derived.  The 
resultant  foim  contains  two  distinct  work  rate  terms  for  the  axial  and  radial  components  of  the  flow.  The  fact  that  integration  yields  a 
result  which  approximates  the  classic  Euler  Turbomachineiy  Equation  constitutes  confirmation  of  the  derivation. 

An  application  of  the  equation  to  an  ideal  infinite  linear  cylinder  with  bound  voiticity  was  developed,  the  cylinder  was  made  to  act 
like  a  turbine  blade  performing  work  by  lifting  an  ideal  airframe  against  gravity.  The  integration  yielded  the  expected  known  result. 


20.  DWTRIBUTION/AVAI LABILITY  OF  AB8TRACT 

SI  UNCLASS IFIEDAJNUMfTED  □  SAME  A8  F 


22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Herman  B.  Urbach 


DO  FORM  1478,  JUN  U 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  users  Unclassified 


22b.  TELEPHONE  (Muda  Ana  Coda) 

301-267-2864 


Prwmus  aOtiens  »pa  obaoleto.  rphi  irtty  ni 


220.  OFFICE  SYMBOL 

Code  2720 


SECURITY  GLASSIFICATION  OF  THIS  PAGE 
UNCLASSIFIED 


01 02-LF-01 4-6602 


CONTENTS 


Page 


Nomenclature  .  iv 

Abstract .  1 

Administrative  Information  .  1 

Acknowledgment  .  1 

Introduction  .  1 

The  Galilean  IVansfonnation  .  2 

The  Interpretation  of  Time  Derivatives  .  4 

Time  Dependence  and  Frame  of  Reference .  4 

Derivation  of  the  Differential  and  Integral  Forms  .  7 

Aspects  of  Time-Dependence  in  the  Absolute  and  Moving  Frames  .  7 

The  Coupled  Substantial  Total  Enthalpy  and  Entropic  Energy  Rate 

in  the  Moving  and  Absolute  Frames .  7 

Uncoupling  the  Substantial  Derivative  of  the  Total  Rothalpy 

and  Entropic  Energy  Rates .  8 

Integration  of  the  Total  Enthalpy  Rate .  10 

Derivation  of  the  Integral  Form  from  the  Differential  Form  .  10 

The  Substantial  Total  Enthalpy  Rate  in  a  Two-Dimensional  Device .  13 

The  Total  Enthalpy  Transfer  Rate  Based  on  Aerodynamics .  13 

The  Stream  Function,  Velocity  and  Relative  Enthalpy  in  the 

Frame  of  the  Blade .  15 

The  Substantial  Total  Enthalpy  Derivative  Without  Rotation  .  15 

References .  19 


FIGURES 

1 .  Configurational  relationships  between  the  absolute  and  moving 

coordinate  systems .  3 

2.  Interpretation  of  the  time-dependent  pressure  term  in  rotor  flow: 

(a)  absolute  frame,  (b)  moving  frame,  and  (c)  graphical  representation 

of  the  pressure  over  time  and  angle  0W. .  5 

3.  The  single-blade  linear  turbine:  (a)  sailplane  version,  (b)  sailboat  version, 

and  (c)  the  rotating  cylinder  blade .  14 

4.  Zone  of  interaction  between  the  fluid  and  the  blade .  16 


iii 


DTRC-PAS-90— 45 


NOMENCLATURE 


/ 

h 

K 

Kjw 

Ho 

L 

m 

n 

P 

Po 

<1 

<1 

r 

n 

s 

t 

T 

U 

V 

(VIS) 

w 

z 

r 


A  unit  vector 
Any  continuous  function 
Enthalpy  per  unit  mass  of  fluid 
Total  enthalpy  per  unit  mass  of  fluid 
Relative  total  specific  rothalpy 

Total  enthalpy  of  a  fluid;  in  linear  systems,  the  total  enthalpy  of  a  fluid 
per  unit  length  of  blade 
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Mass  flow  rate 

A  constant 
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ABSTRACT 


A  differential  turbomachinery  equation  describing  the  energy  transfer  between  a 
fluid  and  any  body  moving  in  that  fluid  was  derived.  The  derivation  is  based  upon  the 
Coriolis  form  of  the  Navier -Stokes  Equations.  A  differential  equation  for  the  total  relative 
rothalpy  was  also  obtained.  The  equalities  contain  a  rigorous  viscous  correction  for  the 
total  enthalpy  and  rothalpy. 

The  differential  equation  defines  the  substantial  derivative  of  the  total  enthalpy  at 
any  point  in  a  fluid.  This  substantial  derivative  represents  the  energy  transfer  rate  into  or 
out  of  the  fluid  at  the  point.  In  order  to  obtain  the  total  enthalpy  transfer  between  the  fluid 
and  the  body,  it  is  necessary  to  integrate  over  the  entire  domain.  For  ideal  flow  regimes, 
the  integration  domain  may  be  restricted  entirely  within  the  rotor.  For  real  fluids,  viscous 
coupling  requires  that  regions  outside  the  blade  area  be  considered. 

On  integration  of  the  differential  equations,  a  form  of  the  Euler  Turbomachinery 
Equation  with  viscous  correction  is  derived.  The  resultant  form  contains  two  distinct  work 
rate  terms  for  the  axial  and  radial  components  of  the  flow.  The  fact  that  integration  yields 
a  result  which  approximates  the  classic  Euler  Turbomachinery  Equation  constitutes  con¬ 
firmation  of  the  derivation. 

An  application  of  the  equation  to  an  ideal  infinite  linear  cylinder  with  bound  vortic- 
ity  was  developed.  The  cylinder  was  made  to  act  like  a  turbine  blade  performing  work  by 
lifting  an  ideal  airframe  against  gravity.  The  integration  yielded  the  expected  known  re¬ 
sult. 
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INTRODUCTION 

A  monument  of  turbomachinery  technology  is  Euler’s  Turbomachinery  Equation 
which  is  based  upon  thermodynamic  definitions  of  work  and  Newton’s  Laws.  Since  the 
Navier-Stokes  Equations  and  Crocco’s  Equation1,2  in  a  rotating  (moving)  frame  are  also 
based  upon  thermodynamics  and  Newton’s  Laws,  they  must  in  principle  contain  Euler’s 
Turbomachinery  Equation  in  differential  form  and,  on  integration,  in  integral  form. 
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Integration  of  Crocco’s  Equation  in  the  absolute  frame  has  given  rise  to  the  “Un¬ 
steadiness  Paradox”  to  explain  rotor  energy  transfer.3"14  An  aspect  of  the  paradox  which 
must  be  rationalized  (vide  infra)  is  Vavra’s  statement  (ref.  8,  pp.  110  and  111)  that  if  the 
relative  velocity  and  the  rotor  velocity  are  constant  with  time,  then  the  absolute  velocity 
must  be  time  independent. 

Time-dependence  is  largely  ignored  in  the  design  of  marine  propellers.  The  fact  that 
energy  transfer  is  routinely  calculated  in  a  steady-state  moving  frame  provides  philosoph¬ 
ical  questions  concerning  the  proper  interpretations  of  the  Unsteadiness  Paradox.  This 
rationalization  is  effected  quite  simply  by  the  “crypto-steady"  relation  that  arises  from  the 
Galilean  transformation  used  by  Coriolis. 

A  Galilean  transformation  with  rotation  (the  Coriolis  transformation)  that  connects 
the  moving  rotor  frame  and  the  absolute  or  laboratory  frame  provides  a  relationship  be¬ 
tween  the  frames  so  that  integration  of  the  energy  rate  may  be  conveniently  performed  in 
a  time-independent  frame  with  a  time-independent  set  of  coordinates. 

The  transformation  leads  to  simplified  expressions  for  the  substantial  total-enthalpy 
transfer  rate  which  are  uncoupled  from  the  expressions  for  the  substantial  entropic  energy 
rate.  The  fact  that  the  total  enthalpy  and  the  entropic  energy  rates  are  uncoupled  makes 
for  simplified  integration  of  the  total  enthalpy  transfer,  and  generates  a  differential  form 
of  the  turbomachinery  equation  and,  on  integration,  a  novel  form  of  Euler’s  Turboma¬ 
chinery  Equation  corrected  for  viscous  non-ideal  flow. 

Euler’s  classical  Turbomachinery  Equation  does  not  include  a  viscous  term  to  cor¬ 
rect  for  viscous  losses  since  Euler  was  concerned  only  with  the  mechanical  work  transfer 
to  or  from  the  shaft.  If  one  is  concerned  with  total  enthalpy  transfer  at  a  point  in  a  fluid, 
the  viscous  loss  must  be  addressed  as  is  illustrated  in  this  paper. 

An  application  of  the  differential  turbomachinery  equation  is  described  for  a  two- 
dimensional,  ideal,  linear  turbine. 

THE  GALILEAN  TRANSFORMATION 

The  moving  frame  and  the  absolute  frame  of  a  turborotor  are  connected  by  a  Gali¬ 
lean  transformation,  which  imposes  relationships  between  the  coordinates  of  the  frames. 
From  these  relationships  the  vector  operations  in  the  two  frames  may  be  derived. 

In  the  following  discussion  the  subscripts  v  and  w  represent  the  absolute  and  the 
moving  frame  coordinate  and  vector  values.  (See  Fig.  1  and  the  Nomenclature  for  defini¬ 
tions  of  quantities.) 

Following  Spannhake,3  derivatives  in  the  moving  coordinate  and  time  have  the  fol¬ 
lowing  relationship: 


Similar  results  are  obtained  with  cartesian  coordinates. 
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The  Galilean  Transformation: 


Fig.  1 .  Configurational  relationships  between  the  absolute  and  moving  coordinate  systems. 


Equations  1  and  2  define,  in  fact,  the  “crypto-steady  criterion.”  15  If  U  is  constant,  then  a 
frame  exists  in  which  the  flow  regime  is  truly  steady  state. 

The  vector  operator  V  is  independent  of  time.  Therefore 

Vv  =  VH,  =  V  .  (3) 

Now  from  Eq.  1  for  any  function/, 
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where  4  represents  all  the  position  coordinates  in  the  absolute  frame  and  qwu  is  the  mov¬ 
ing-frame  coordinate  of  a  generalized  curvilinear  orthogonal  coordinate  system  which 
completely  defines  the  position  and  velocity  of  the  energy-transferring  device.  Eq.  4  is  an 
extension  of  the  usual  crypto-steady  relation. 

THE  INTERPRETATION  OF  TIME  DERIVATIVES 

The  time  derivative  of  the  static  blade-to-blade  pressure  is  obtained  from  Eq.  4  by 
substituting  p  for  the  function /.  Using  cylindrical  coordinates 


Figure  2a  shows  a  point  P(r,z,t )  fixed  in  the  absolute  frame  between  blades  of  a  cen¬ 
trifugal  compressor  rotor  rotating  into  decreasing  values  of  6V  .  The  thermodynamic 
properties  of  the  point  P  at  the  suction  side  of  blade  1  change  as  time  advances  and  the 
pressure  side  of  blade  2  approaches  P.  When  the  wall  passes  through  point  P,  fluid  prop¬ 
erties  cease  to  exist  at  P.  Thermodynamic  information  about  the  fluid  at  this  point  ceases. 
The  time  interval  At  for  n  blades  and  angular  velocity  (1)  is 

t2-t\-  Mr, z)  =  [2 Ji/n -  <5(r , z)]/ty  ,  (6) 

where  <5(r,  z)  is  the  subtended  blade  angle  in  radians. 

Tune  is  fixed  at  to  for  the  spatial  derivatives  of  Fig.  2b.  The  angular  change  A Qw  is 
in  a  direction  negative  to  the  blade  motion. 

Figure  2c  shows  a  curve  of  the  intra-blade  pressure  distribution  changing  with  time 
in  the  absolute  frame.  In  the  moving  frame,  it  represents  the  spatial  distribution  in  a  fixed 
instant  of  time. 


aA  =_a(?p\ 

3,)r^  r\S0»)W 


dp 


d(6w/(o)t 


(5) 


TIME  DEPENDENCE  AND  FRAME  OF  REFERENCE 

Looking  at  the  relationship  between  moving  and  absolute  frame  velocities,  namely 

V=W+U ,  (7) 

Vavra  has  noted  that  Eq.  7  suggests  that  if  W  is  independent  of  time,  then  V  is  also  inde¬ 
pendent  of  time.  This  point  is  an  overlooked  aspect  of  the  Unsteadiness  Paradox.3-14  The 
converse  statement  is  also  true,  and  the  Unsteadiness  Paradox  would  imply  that  flow  must 
be  unsteady  in  the  moving  frame.  The  observer  who  sits  on  an  ideal  rotor  in  an  ideal  infi¬ 
nite  fluid  sees  no  time  dependence  in  measured  thermodynamic  properties  at  a  point. 
However,  when  the  same  observer  passes  to  the  absolute  frame,  he  measures  time-depen¬ 
dent  thermodynamic  properties  each  time  a  blade  passes  by  that  point. 

The  problem  is  resolved  by  noting  that  the  frame  of  the  coordinates  used  in  Eq.  7 
determines  whether  the  observer  perceives  time  dependence.  If  the  velocity  vectors  are 

expressed  in  terms  of  r,  6W,  and  z,  the  measurements  are  in  the  moving  frame  and  both 

V(r,  6W,  z)  and  W(r,  6W,  z)  are  time  independent.  Writing  V  and  W  in  the  absolute 
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a)  Absolute  frame 

P(r,  0V,  z)  is  fixed  in  space.  Time  changes.  Blades  rotate. 

SUCTION  SIDE 


£p\  =  —  (  ep  \ 

/  r,0v.  z  r  \  e0w  /  r,  z,  t 


OVER  TIME,  M  THE  PRESSURE  SIDE  OF  THE  NEXT  BLADE 
MOVES  TO  THE  POINT  P(r,  0v,  z). 


Fig.  2.  Interpretation  of  the  time-dependent  pressure  term  in  rotor  flow:  (a)  absolute  frame,  (b)  moving 
frame,  and  (c)  graphical  representation  of  the  pressure  over  time  and  angle  6W. 
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PRESSURE 


b)  Moving  frame 

P(r,  2,  to)  over  a  range  of  decreasing  0W  in  a  fixed  instant,  t0.  Blades  appear  stationary. 


c)  Graphical  representation  of  the  pressure  over  time  and  angle  0W. 
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coordinates  requires  that  6W  be  written  6v-(ut  .  Therefore,  V(r,6v-a)t,z)  and 
W(r, 6v-(Ot ,  z)  are  both  time  dependent  in  the  absolute  frame. 


DERIVATION  OF  THE  DIFFERENTIAL  AND  INTEGRAL  FORMS 


ASPECTS  OF  TIME-DEPENDENCE  IN  THE  ABSOLUTE  AND  MOVING  FRAMES 

The  Navier-Stokes  Equations  exhibit  Coriolis  terms  in  a  unique  steady-state  rotating 
frame  which  develops  from  the  Galilean  transformation.  Even  in  the  frames  of  a  device 
without  rotation,  the  definition  of  the  moving  frame  inherently  contains  and  always  in¬ 
vokes  the  crypto-steady  relation,  Eqs.  4  and  5.  Provided  there  is  uniform  upstream  flow, 
an  observer  moving  with  the  device  (for  example  a  sail)  may  observe  crypto-steady-state 
behavior  in  the  fluid.  Nevertheless,  the  time-dependent  terms  must  be  recovered  on  impo¬ 
sition  of  the  Galilean  transformation  to  the  absolute  frame. 

In  the  Coriolis  form  the  Navier-Stokes  Equations8  explicitly  exhibit  the  rotational 
motion  of  the  energy-transferring  rotor  device,  thus: 

3(7  _  _  ^  $7  _  _ 

—  +V-VV  = - +  WVW  +  2a>x  W-VU2/ 2 

dt  dt  ' 


Vp  1  - 

=  — — +  — V-jt 


(8) 


where  jr'  represents  the  stress  tensor  excluding  the  pressure  term  pd 


j 
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THE  COUPLED  SUBSTANTIAL  TOTAL  ENTHALPY  AND  ENTROPIC  ENERGY 
RATE  IN  THE  MOVING  AND  ABSOLUTE  FRAMES 

If  crypto-steady  flow  does  not  characterize  the  flow  regime  of  the  moving  frame, 
the  thermodynamic  properties  of  the  flow  must  fluctuate  about  an  average  value  for  any 

given  thermodynamic  condition.  If  U  is  constant,  the  time-dependent,  linear  acceleration 
term  evaluated  in  the  relative  frame  integrated  over  time  cannot  lead  to  energy  transfer. 
However,  time  averages  over  non-linear  terms  lead  to  non-vanishing  Reynold’s  stress 
terms.  It  will  be  assumed  that  time  dependence  in  the  moving  frame  is  negligible.  There¬ 
fore  in  Eq.  8  the  acceleration  term  in  the  moving  frame  will  be  dropped. 

~  +  W2/2-V  x  (V  x  V) 


=  VW2p  -W  x  (V  x  W)  +  Wx  W-VU2/ 2 


-  Vp  1  — 

- -  +  —  V -7t' 


(9) 


Using  the  absolute-frame  equality  of  Eq.  9  and  the  gradient  form  of  Gibb’s  equation 
of  state,  the  substantial  derivative  of  the  total  enthalpy  is  obtained  in  terms  of  the  partial 
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Ds 

derivative  of  the  pressure,  the  substantial  entropic  energy,  T ,  and  the  stress  tensor 


7t'  (see  refs.  16-18  for  deta'ls.) 

Dh0  1 
Dt  Q 


+  r— +-  (V-P) . 

Dt  Q 


(10) 


If  the  substantial  entropic  energy  is  independent  of  heat  transfer,  then  it  may  be  ex¬ 
pressed  as  functions  of  the  velocity  vector  along  with  the  stress  term.  At  this  point 
analysis  stops  unless  the  partial  pressure  derivative  is  resolved  in  terms  of  the  velocity. 
Employing  Eq.  4  with  the  static  pressure  as  the  arbitrary  function  the  result  is 

Dh0  V  Ds  V 

— - -  Vn  +  7 - 1 - -  i  v  •  jt  i 

(ID 


—  - - ■Vp  +  T~^  +  ~  •  (V ■  jt’)  . 

Dt  Q  Dt  g 


where  Eq.  9  may  be  used  to  replace  the  pressure  derivative  in  Eq.  11.  From  the  dot  prod¬ 
uct  of  U  with  the  moving  frame  equality  in  Eq.  9  the  pressure  gradient  term  becomes 

—  =  -£/•  —  =  Z7  •  (W  •  V)VF  +  2(7  •  ST  X  W 

Qr  d0„  Q 


-  U  ■  VU2/2  — —  1/  •  (V  •  n') 
Q 


(12) 


In  Eq.  12  the  gradient  of  the  pressure  is  of  course  invariant  in  all  frames,  but  the 
moving  frame  is  convenient  and  preferable.  Combining  Eq.  11  with  12,  the  substantial 
derivative  of  the  total  enthalpy  coupled  with  the  entropic  energy  is  obtained  in  terms  of 
flow  field  variables. 


Dh, 


.  Ds  „ 


1 


~2.-T—=U  (W  V)W  +  2U  a)X  w  +  -W-(y  -n')  . 
Dt  Dt  g 


(13) 


UNCOUPLING  THE  SUBSTANTIAL  DERT.ATTVE  OF  THE  TOTAL  ROTHALPY 
AND  ENTROPIC  ENERGY  RATES 

If  the  gradient  of  the  specific  kinetic  energy  in  the  absolute  frame  is  subtracted  from 
the  right-hand  equality  in  Eq.  9  and  V  is  replaced  by  Eq.  7,  then18 

-VWUU-VU2  +  (V  x  W  +  2g7)  x  W=-Vh0  +  TVs+-V -P  . 

Q  (14) 


A  useful  intermediate  expression  is 

v2  ,  w2  u2 

=  h  +  —  -UWu-u2  =  h  +  —  --j  , 


(15) 


where  the  quantity  h ^  in  Eq.  15  is  the  total  relative  rothalpy.  Combining  Eqs.  14  and  15, 
Crocco’s  Equation  for  the  total  relative  rothalpy  in  the  moving  frame  with  crypto-steady 
flow  is  derived. 
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(16) 


Vhow-TVs  =  W  x  (V  x  W  +  2fiJ)  +— Vj?  . 


Now,  since  the  upstream  flow  is  thoroughly  mixed  and  without  energy  gradients,  time 
derivatives  of  the  pressure,  rothalpy  and  entropy  vanish  in  the  moving  frame.  Thus 

Dt  dt  (17) 


Jds\  =J,d1 

\*)w  * 


+  TW-Vs=lW'Vs  =  T 


where  the  last  equality  in  Eq.  18  indicates  that  the  entropic  energy  rate  is  invariant  in  all 
frames.  Ihking  the  dot  product  of  W  on  Eq.  16  yields 

—  =  r^+-i^-(V  P) . 

Dt  Dt  g  (19) 

Now  Wu  notes  that  in  the  moving  frame  following  a  particle  of  fluid  (ref.  17,  p.  91 ) 

rfHV*'V)^  (20) 

where  q  arises  from  heat  transfer  and  the  second  term  is  the  specific  dissipation.  It  is  as¬ 
sumed  that  there  are  no  external  heat  sources  or  sinks  and  no  thermal  conduction  and  that 
crypto-steady  flow  prevails.  Therefore 

T^  =  -(?  V)  W  =  -<p  , 

Dt  q  g  (21) 

where  <f)  is  the  dissipation.  Combining  Eqs.  19  through  21,  the  desired  result  is 

Dt  Q  (22) 

Equation  22  asserts  that  the  substantial  derivative  of  the  total  relative  rothalpy  is 
dependent  upon  the  viscous  dissipation.  Using  Eqs.  17  and  22  for  inviscid  flows  without 
heat  sources  or  sinks,  and  with  crypto-steady  characteristics 


V/inR,  =  V I  h 


W2  U2)  _ 
2  ~  2  )  ~ 


V(ho-UVu)  =  0  . 


Equation  23  represents  differential  or  gradient  forms  of  Euler’s  Turbomachinery  Equa¬ 
tions  in  the  moving  and  absolute  frame.  Integration  over  a  stream  tube  in  the  absolute 
frame  yields  the  classic  Euler  Equation. 
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A h0  =  A(UVU)  . 


(24) 


Note  that  Eq.  24  is  strictly  true  only  for  isentropic  flow.  The  substantial  total  enthalpy  rate 
in  the  absolute  frame  is  now  obtained  by  eliminating  the  entropic  energy  rate  in  Eq.  13 
with  Eq.  21. 

=  U- (W •  V)fV  +  2U ■  (O  x  W  +  - V  •  (n'  ■  W)  . 

Dt  Q  (25) 

Equation  25,  expressing  the  power  transfer,  is  a  differential  form  of  the  turbomachinery 
equation  in  terms  of  the  moving  frame. 

With  Eq.  21  it  is  now  possible  to  uncouple  the  total-enthalpy,  pressure  relationship 
of  (11)  from  the  entropic  energy. 

^  =  _^.Vp  +  iv-(^-W0+if7-(V ■*')  . 

Dt  e  Q  Q  (26) 

Equation  26  indicates  that  neither  axial  nor  radial  pressure  gradients  are  germane  to 
the  calculation  of  specific  total  enthalpy  transfer.  (Mass  flow  rates  are  of  course  a  func¬ 
tion  of  axial  or  radial  pressure  gradients.)  Note,  in  contrast  with  some  views  (ref  9,  pp.  7 

and  8),  only  transverse  pressure  gradients  parallel  to  U  contribute  to  total  enthalpy  trans¬ 
fer.  The  impulse  stages  of  turbomachines  prove  that  axial  or  radial  pressure  gradients  play 
no  role  in  energy  transfer.  Moreover,  from  Eq.  25  for  ideal  flows  which  have  no  vorticity, 

only  kinetic  energy  gradients  paralleling  the  blade  motion  U  contribute  to  total  enthalpy 
transfer.  Nevertheless,  non-ideal  total  enthalpy  transfer  depends  upon  non-linear  as  well 
as  the  linear  properties  of  the  flow. 

The  substantial  total  enthalpy  rate  given  by  Eq.  25  is  a  differential  form  of  the  tur¬ 
bomachinery  equation  in  the  time-independent  coordinates  of  the  moving  frame.  A 
proper  test  of  Eq.  25  would  be  the  applicability  of  the  equation  to  integration  over  the 
rotor  blade-to-blade  flow.  Moving-frame  integration  should  predict  a  total  enthalpy  trans¬ 
fer  compatible  with  that  of  the  Euler  Turbomachinery  Equation.  Therefore,  the 
integration  of  Eq.  25  will  be  performed  as  a  test  in  the  three-dimensional  domain.  An 
application  of  the  new  equation  will  be  developed  in  a  two-dimensional  linear  turbine. 
Compatibility  of  the  results  with  Euler’s  Equation  or  the  consequences  thereof  will  lend 
credence  to  the  logic  of  analysis  employed  in  the  derivation. 

INTEGRATION  OF  THE  TOTAL  ENTHALPY  RATE 

In  the  integration  process  it  will  be  assumed  that  the  flow  may  be  divided  into 
streams  which  pass  between  a  given  pair  of  blades.  In  the  rotating  frame  the  stream  tube 
walls  are  fixed  steady-stage  walls  associated  with  a  steady-state  mass  flow  rate  m  which 
may  consist  of  radial  and  axial  mass  flow  components. 

DERIVATION  OF  THE  INTEGRAL  FORM  FROM  THE  DIFFERENTIAL  FORM 
The  differential  form  of  the  turbomachinery  equation  (Eq.  25)  may  be  integrated  to 

yield 
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Q  [U  ■  (W  •  V)W+2U  -ioxw]  rdrdQdz 


V  •  (j?  •  wo  rdrdfldz  . 


The  first  term  of  the  right  member  of  Eq.  27  is  the  tangential  component  of  the  convec¬ 
tive  term  obtained  on  dot  multiplication  with  U  ,  i.e.: 

eU-(W-  V)W 


mtUlWrZt+*«*L.w.Z!!-+Z&). 

^  \  dr  r  86  dz  r  ) 


The  first  and  last  terms  of  Eq.  28  will  be  combined  in  an  integral  identified  by  / 1,4  thus: 

hA  =\\\erW~i  ^-rdrdddz  .  (2( 

The  factors  in  Eq.  29  may  be  rearranged  to  express  the  radial  mass  flow  which  must  be 
constant  in  steady  flow. 

Following  arguments  presented  previously18 

/1>4  =  mr  J  d(UW0 )  =  mr[{lf^g)r2  -  (U&eU]  ,  (3( 

where 


The  third  term  of  the  right  member  of  Eq.  28  may  be  written  to  show  the  axial  mass 
flow  rate  m*  explicitly. 

h  =  /(/ / e^*rdrde  j ^~~dz  •  (32) 

where  the  axial  velocity  Vz ,  Wq  ,  and  U  have  been  averaged  over  A r  and  Ad  .  Substi¬ 
tuting  the  axial  mass  rate  for  the  parenthesis  in  Eq.  32, 

h  =  mz  j  d(UWe )  =  mz[(UWe)z 2  -  (0W)Z  1]  .  (33) 
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The  second  term  of  the  right  member  provides  an  integral  which  contains  the  tan¬ 
gential  kinetic  energy. 


30  / 


=  j  - WfiO^drdz  =  0 


Since  the  tangential  velocities  at  the  blade  walls  are  equal  to  the  blade  velocity,  the 
integral  vanishes. 

Now  identifying  the  second  term  of  the  right  member  of  Eq.  27  as  I5,  we  may  write 


•  . 


2Q(D2rWrrdrd6dz  . 


Finally,  following  the  arguments  above, 

^5  «  Wr[(l/2)r2-(tf2)rl]  • 

Summing  the  components  of  integration  /jf4  through  1$, 


dk.  rA 

1 - dr  =  mr  Ari 

Dt 


(UWe)  +  AA.U2) 


+  mzAz(UW0)  +  (VIS) 

where  Ar  and  Az  represent  the  change  along  r  and  z  respectively,  and  (VIS)  is  the  inte¬ 
grated  viscous  term. 

Now  adding  m2Az(72,  which  is  zero,  to  Fq.  37, 


JJJe^rfr=mrAr  (U+W0)U 


+  m2 Az  (U  +  We)U\+(VIS)  . 


In  Eq.  38  the  terms  Wg  are  averaged  over  6  and  z  in  the  Fust  term  and  over  r  and 
6  in  the  second  term.  If  the  total  steady-state  mass  rate  m  between  a  pair  of  blades  is 

m  =  mr  +  mz  ,  (39; 
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A ho  =frAr(UVe)  +f:Az(UVe)  +  (VIS)/m  .  (40) 

Equation  40  represents  Euler’s  Turbomachinery  Equation  with  mixed  flows  and  the 

terms  Vg  and  U  are  averaged  over  the  blade  space  where  necessary.  The  coefficients  fr 
and  f  represent  the  radial  and  axial  fractions  of  the  mass  flow. 

The  viscous  term  is  a  novel  feature  of  the  derivation  which  may  explain  losses  of 
rotors  during  windmilling.  With  the  exception  of  the  viscous  energy  term  (VIS),  the  inte¬ 
gral  expression  Eq.  40  exhibits  a  formal  similarity  and  compatibility  with  Euler’s 
Turbomachinery  Equation.  The  derivation  lends  credence  to  the  hypothesis  that  Eq.  25  is 
indeed  a  differential  form  of  the  turbomachinery  equation  (Eq.  24). 

A  two-dimensional  application  and  test  of  the  differential  form  (Eq.  25)  on  an  ideal 
linear  device  where  the  solution  is  known  precisely  will  now  be  examined. 

THE  SUBSTANTIAL  TOTAL  ENTHALPY  RATE  IN  A  TWO-DIMENSIONAL 

DEVICE 

An  infinite  circular  cylinder  with  bound  circulation,  as  shown  in  Fig.  3,  is  an  ele¬ 
mental  linear  turbine.  It  may  be  considered  as  an  infinite  sail  on  a  sailboat  or  an  infinite 
wing  on  a  sailplane.  The  device  extracts  energy  from  the  ideal  inviscid  working  fluid. 
Work  is  performed  on  the  sailplane  (fixed  to  a  vertical  rail)  by  raising  its  height  at  uni¬ 
form  speed  U  against  gravity.  Work  on  the  sailboat  is  performed  by  moving  the  boat  at 
uniform  speed  U  which  elevates  a  weight  attached  at  minus  infinity  by  an  infinite  tether. 

In  the  moving  frame  the  apparent  velocity  of  the  ideal  working  fluid  at  infinite  distance  is 
W0  .  The  relationship  between  the  absolute  and  moving  coordinate  system  and  the  velo¬ 
cities  is  given  by  the  Galilean  transformation  of  Fig.  1 . 

THE  TOTAL  ENTHALPY  TRANSFER  RATE  BASED  ON  AERODYNAMICS 

In  the  absence  of  vorticity  and  viscosity  the  ideal  flow  field  will  exhibit  no  total  gra¬ 
dient  in  accordance  with  Eq.  23.  Since  the  flow  field  is  ideal,  the  flow  domain  may  be 
described  by  a  potential  function  or  its  conjugate  stream  function.  The  lift  is  therefore  the 
ideal  lifting  force,  L,  of  the  Kutta-Joukowski  Equation  given  by 

L  =  eW<£  ,  (41) 

where  T  is  the  scalar  circulation.  The  units  are  force  per  unit  length  of  cylinder.  In  the 
absolute  and  moving  frame  the  lift  component  Ly  directed  parallel  to  the  y  axis  of  Fig.  3  is 
given  by 

Ly  —  qW0XT  —  qVoxT  ,  (42) 

where  the  subscript  x  represents  the  x  component.  Recalling  that  U  is  the  velocity  of  mo¬ 
tion  of  the  device  (sail  or  wing  or  rotating  cylinder)  as  perceived  in  the  absolute  frame, 
the  power  is  the  product  of  U  and  Ly. 

Power /unit  length  =  qUW0XT  .  (43) 

Since  we  assume  that  there  is  no  heat  rate, 


DTRC-PAS-90-45 


13 


(44) 


DH0 

Dt 


=  -qUWoxT 


where  H0  is  the  total  enthalpy  of  the  system  per  unit  length.  Equation  44  is  the  anticipated 
relationship  which  should  ultimately  be  developed  from  the  differential  form  Eq.  25. 


THE  STREAM  FUNCTION,  VELOCITY  AND  RELATIVE 
ENTHALPY  IN  THE  FRAME  OF  THE  BLADE 

Since  ideal  flow  has  been  assumed  in  the  moving  frame  of  the  blade,  the  stream 
function,  \p  ,  is  the  usual  function  modified  for  motion  along  the  y  axis. 

xp  =  -  Woy(  1  -  a2/r2)r  cos  6  +  Wox(  1  -  a2/r2)r  sin  6 


+  (T/2jt)ln(r/a)  .  (45) 

The  constant  a  is  the  radius  of  the  cylinder.  The  cartesian  velocity  components  are 
obtained  by  the  usual  transformation18  as  follows: 

w  -w  i ftygjyj 

x~  °x  Or2  +  y2)2  (x2  +  y2)2 

+  _L_Z_ 

2ji  (x2  +  y2)  (46) 

2 a2Woxxy  .  Jw^-y1) 
y~  (x2^/)2  (x2  +  y2)2 

Y_ _ x 

ht  (x2  +  y2)  '  (47) 

Now  the  relative  vorticity  must  vanish  because  potential  flow  cannot  have  vorticity. 
A  check  of  the  vorticity  in  the  relative  frame  shows  that  indeed  it  vanishes.  Also,  the 
time-dependent  term  vanishes. 


THE  SUBSTANTIAL  TOTAL  ENTHALPY  DERIVATIVE  WITHOUT  ROTATION 


In  the  linear  two-dimensional  system,  the  differential  form  of  the  turbomachinery 
equation  (Eq.  25)  is  simplified  because  the  rotation  vanishes. 

~-=U  (W  V)W  .  (48) 


Since  the  vorticity  vanishes 


W-VW=VW2/ 2  , 


and  employing  Eq.  4 


^=(7.V^/2  =  -—  • 

Dt  '  2  dyw 


DTRG-PAS-90-45 


15 


Fig.  4.  Zone  of  interaction  between  the  fluid  and  the  blade. 


The  integrated  substantial  total  enthalpy  rate  per  unit  length  (where  the  subscript  on  y  in 
Eq.  50  has  been  dropped)  is 


BUz.  I  [  «£  wUvfy 

Dt  J  J  2  dy  3 


(51) 


Integration  of  Eq.  51  will  be  performed  over  all  space  per  unit  length  z  of  the  blade.  The 
choice  of  time  is  immaterial  since  the  fluid  dynamics  are  steady  state  in  the  moving  frame 
and  the  thermodynamic  rates  over  all  space  are  invariant  with  time.  It  is  understood  that 
the  integration  applies  only  to  the  fluid  domain  and  that  boundaries  at  solid  walls  are  ob¬ 
served. 


r  00  r  so 

=  qU  I  d(W2/l)dx  . 

J  -00  J  -00 


See  Fig.  4  for  a  definition  of  the  boundary  conditions  of  y  on  the  blade  wall.  Since  the 
kinetic  energy  in  a  conservative  system,  like  thermodynamic  quantities,  is  a  function  of 
state,  the  integration  in  Eq.  52  is  a  function  of  the  end  points  only,  which  is  verified  by 
the  fact  that  the  integrand  is  a  total  differential. 
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Since  the  velocity  is  uniform  at  infinity  (and  all  time  and  space  derivatives  vanish  at  in¬ 
finity),  the  contributions  at  infinity  will  cancel  in  Eq.  53. 


According  to  Fig.  4  integration  of  Eq.  53  with  respect  to  x  yields  contributions 
from- <,  x  :£  a  only. 


DH0 

Dt 


ra 

=  -qV\  d(Wl  +  Wy)/2dx  . 

J -a  J -(a1-*2)''2 


(54) 


Now  without  going  into  details18  the  integration  of  the  total  derivative  in  Eq.  54  is  per¬ 
formed,  yielding 


DH0 

Dt 


~  —qU 


4 UWox  Jtfa2-*2)1/2^-*2) 


IWpxTxtf-x2)1/2  &WoxUx(a2  -  x2)1^2 
7ta4  a2 


WoxT(a2  -  jc2)1/2  AW0XUx{a2  - x 2)ll\a2  -  lx2) 

+  ^  + - 7 - 

Woxr(a2-x2)l/2(a2-lx2)\j 

+ - ^ (55) 

Note  that  only  odd  terms  in  y  make  any  contribution  to  Eq.  55.  Since  the  first  and  fifth 
terms  cancel,  only  four  terms  remain.  The  integration  with  respect  to  x  is  performed 
through  a  transformation  employing 

x  =  acos6  ,  (56) 

with  integration  limits  given  by 

6  =  71  when  x  =  -  a  , 

0  =  0  when  x-  a  .  (S'7) 


Making  the  substitutions 
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(58) 


In  Eq.  58  the  second  integral  makes  no  contribution  because  it  is  antisymmetric.  The  first 
integral  cancels  the  second  term  in  the  last  integral  to  yield  from  the  surviving  terms 

DH 

-  QWoxrU  =  -  QVoxTU  .  (59) 

Equation  59  is  identical  with  Eq.  44  and  this  result  illustrates  a  useful  application  of 
the  differential  form  and  constitutes  confirmation  of  the  validity  of  the  differential  turbo¬ 
machinery  Equation  (Eq.  25).  For  the  linear  case,  the  energy  transfer  rate  of  the  rotor  is 
proportional  to  the  component  of  the  kinetic  energy  gradient  parallel  to  the  moving  rotor 
(or  sail). 

Recall  that  energy  transfer  occurs  in  the  narrow  domain  of  integration  indicated  in 
Fig.  4  which  ranges  to  infinity.  Thus  the  velocity  of  sound  must  be  infinite  in  the  potential 
system  in  agreement  with  the  assumption  of  incompressibility. 
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